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Abstract 

Background: Research on the neural bases of cognitive deficits in autism spectrum disorder (ASD) has shown that 
working memory (WM) difficulties are associated with abnormalities in the prefrontal cortex. However, cognitive 
load impacts these findings, and no studies have examined the relation between WM load and neural 
underpinnings in children with ASD. Thus, the current study determined the effects of cognitive load on WM, using 
a visuo-spatial WM capacity task in children with and without ASD with functional magnetic resonance imaging 
(fMRI). 

Methods: We used fMRI and a 1-back colour matching task (CMT) task with four levels of difficulty to compare the 
cortical activation patterns associated with WM in children (7-13 years old) with high functioning autism (N= 19) 
and matched controls [N= 17) across cognitive load. 

Results: Performance on CMT was comparable between groups, with the exception of one difficulty level. Using 
linear trend analyses, the control group showed increasing activation as a function of difficulty level in frontal and 
parietal lobes, particularly between the highest difficulty levels, and decreasing activation as a function of difficulty 
level in the posterior cingulate and medial frontal gyri. In contrast, children with ASD showed increasing activation 
only in posterior brain regions and decreasing activation in the posterior cingulate and medial frontal gyri, as a 
function of difficulty level. Significant differences were found in the precuneus, dorsolateral prefrontal cortex and 
medial premotor cortex, where control children showed greater positive linear relations between cortical activity 
and task difficulty level, particularly at the highest difficulty levels, but children with ASD did not show these trends. 

Conclusions: Children with ASD showed differences in activation in the frontal and parietal lobes — both critical 
substrates for visuo-spatial WM. Our data suggest that children with ASD rely mainly on posterior brain regions 
associated with visual and lower level processing, whereas controls showed activity in frontal lobes related to the 
classic WM network. Findings will help guide future work by localizing areas of vulnerability to developmental 
disturbances. 
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Background 

Approximately 1 in 88 children have an autism spectrum 
disorder (ASD), and the proportion of children being diag- 
nosed is rising [1]. ASD is a neurodevelopmental disorder 
classically characterized by social deficits, communicative 
difficulties and repetitive behaviours [2], with evidence of 
cognitive and executive function impairment [2-9]. The 
impaired executive processes may account for many of the 
profound behavioural manifestations in ASD and contrib- 
ute to autistic symptomology [6,7]. Emerging literature on 
cognitive difficulties in ASD, as well as their neural under- 
pinnings, provides evidence for working memory (WM) 
deficits that are associated with frontal lobe abnormalities, 
particularly in prefrontal cortical activity [4,10-14]. 

Previous neuroimaging, electrophysiology and neuro- 
chemical studies in ASD have identified atypical white 
and grey matter volumes [15-18], functional connectivity 
[10,19], cortical sulcal and gyral anatomy [20], brain 
lateralization [21], neural perfusion [22] and serotonin 
synthesis capacity [23] compared to those of typically 
developing (TD) individuals, with the frontal cortex impli- 
cated in a number of these differences [24]. Although 
there is substantial evidence for developmental anatomical 
abnormalities of the frontal lobes in ASD, associations 
with cognitive performance are lacking [25]. It remains to 
be determined if executive dysfunction in ASD is related 
to functional as well as anatomical frontal lobe aberra- 
tions. Further, protracted maturation of the frontal lobes 
makes the functions they support, such as WM, suscep- 
tible to developmental disturbances [26,27] but amenable 
to therapeutic interventions. 

WM allows for the temporary storage, rehearsal and 
maintenance of information. It is important for learning, 
social ability [28], academic achievement [29] and many 
complex cognitive operations [30,31]. Impairments in WM 
have been reported in ASD, yet our understanding of its 
development and neural correlates is still very limited in 
this population. We investigated the neural systems 
underlying visuo-spatial WM capacity using functional 
magnetic resonance imaging (fMRI) to determine the 
role of cognitive load on possible functional differences 
between children with and without ASD. Cognitive load 
is a multidimensional construct referring to the pro- 
cessing resources that performing a particular task im- 
poses on one's cognitive system [32] — in the present 
study, the level of WM task difficulty was used as a 
measure of cognitive load. 

Extant research on the behavioural characterization of 
WM function in ASD suggests both intact [33-37] WM 
performance on simple memory tasks, and impaired 
[33-36,38] performance on more complex tasks, relative 
to TD individuals. Discrepancies in the literature may be 
due partly to methodological inconsistencies regarding 
task choice, comparison groups or population age. Basic 



WM abilities are intact in high-functioning ASD (see 
[39] for review). The majority of studies that found sig- 
nificant group differences assessed WM using tasks with 
increased complexity and/or cognitive load and thereby 
imposed heavier demands on WM and executive func- 
tions [9,34,38]. Further, there is evidence of primarily 
visuo-spatial WM impairment, whereas verbal WM ap- 
pears relatively intact in individuals with ASD [35-37]. 

Only a few studies have examined the neural correlates 
of WM function in ASD, and no neuroimaging studies 
of pre-adolescent children exist. Overall, WM processes 
are largely subserved by the prefrontal and parietal 
cortices [4,40-44]. Converging literature identifies a 
broad system of prefrontal, premotor, dorsal cingulate 
and posterior parietal activation in visual WM tasks 
see [44], and neuroimaging studies in ASD provide evi- 
dence for atypical activity in these regions [10,12,14]. 
Using an oculomotor visuo-spatial WM task, Luna and 
colleagues [12] found behavioural impairments in 
WM, as well as reduced activation in the dorsolateral 
prefrontal cortex (DLPFC) and posterior cingulate cortex 
in adults with ASD compared to controls. Interestingly, 
there was no evidence of impaired activation in other 
areas known to support WM neural circuitry, such as the 
anterior cingulate cortex. A study by Koshino et al. [10] 
demonstrated bilateral activation in the DLPFC in adult 
controls, whereas adults with ASD showed limited activa- 
tion in the left, and more right hemisphere recruitment of 
prefrontal regions, despite the absence of differences in 
performance accuracy on an n-back letter task. The au- 
thors suggested that while typical adults processed letter 
stimuli using verbal codes, those with ASD employed 
visual strategies [10], supporting the idea that WM defi- 
cits in ASD may be attributed to less efficient process- 
ing strategies. In the only study of adolescents, Silk 
et al. [14] used a visuo-spatial WM mental rotation task, 
where behavioural performance was similar for ASD 
and control groups, and observed impaired cortical acti- 
vation in the frontal lobes in the ASD group, including 
the anterior cingulate, DLPFC and caudate nucleus, but 
normal activation in the parietal cortices relative to con- 
trols. These findings suggest dysfunctional frontostriatal 
networks in ASD. More work is needed to understand 
the neural correlates and developmental trajectory of 
WM in childhood in ASD, as limited research is avail- 
able reporting on pre-adolescent children. 

One of the most common experimental paradigms used 
to manipulate cognitive load in studying WM is the 
'n-back' protocol [10,41,43-47]. In a typical n-back 
task, participants view a series of stimuli and indicate 
whether the currently presented stimulus matches one 
presented 'n' (e.g., 0, 1, 2 or 3) trials prior. As difficulty 
level increases, the number of interfering stimuli between 
the target and relevant stimulus increases, requiring the 
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utilization of different mental strategies at each level (e.g. 
O-back, recognition; 1-back, maintenance; 2-back, main- 
tenance and monitoring). This manipulation consequently 
increases both memory load and executive function 
demand (i.e. strategy needed to solve the task) in a 
non-linear fashion from one level to the next, making 
function-specific changes difficult to quantify and link 
with brain areas. Therefore, we used a 1-back colour 
matching task (CMT) [48] which systematically manip- 
ulated memory load while keeping executive function 
constant across all difficulty levels, allowing a direct 
investigation of the influence of cognitive load on WM. 
Specifically, executive schemes (i.e. procedural strat- 
egies for solving the task) are constant across levels of 
items in CMT; what varies with each level is the number 
of relevant items (colours) to be remembered. Difficulty 
level was parametrically graded based on behavioural age- 
dependent growth patterns observed using this task in 
previous work [48]. These observations point to a linear 
pattern of WM development, which may also be evident 
in the neural processes across increasing load. Our task is 
novel, and it captures the brain regions associated with 
this linear pattern of function across cognitive load in pre- 
adolescent children. In a previous study from our group, 
typical adults showed positive linear relations between 
cortical activity and CMT task difficulty level in areas in- 
volved in WM function [49]. Negative linear relations 
were found in areas typically associated with the default 
mode network (DMN). The DMN, which has been found 
in a wide range of neuroimaging studies, is a network of 
brain regions, including the medial prefrontal cortex, pos- 
terior cingulate and inferior parietal lobules, characterized 
by decreased activation during goal-oriented or attention- 
demanding tasks [50]. Neurodevelopmental disorders, 
including ASD, have been associated with abnormal func- 
tion [51,52] and structure [53] of the DMN that may inter- 
fere with cognitive function. 

The present study identified and compared the neural 
activity underlying visuo-spatial WM capacity in children 
with and without ASD using fMRI, with manipulation of 
cognitive load. Given the increasing difficulty levels of the 
current task, we hypothesized that children with ASD 
would show poorer performance than matched TD con- 
trols on behavioural measures of WM. Further, frontal and 
posterior parietal cortical areas related to visuo-spatial WM 
capacity would be under-recruited in children with ASD 
relative to TD controls, and this difference would increase 
with load. We expected activity to be linearly modulated 
(positively for WM areas, negatively in default-mode areas) 
by task difficulty in all children (comparable to previous 
work with adults), but these trends would be less marked 
in children with ASD, particularly at higher cognitive loads. 

With the protracted maturation of the frontal lobes, 
and their susceptibility to developmental anomalies, 



understanding development in these regions is crucial, 
particularly in populations with frontal lobe abnormal- 
ities, as seen in ASD. As previous research has focused 
on studying WM in adults and adolescents with ASD, 
there is a critical need to study WM in pre-adolescents 
with ASD, who exhibit both cognitive and neurological 
differences. Given the links between social function, 
school success and executive function ability, investigat- 
ing the neural bases of WM deficits in children with 
ASD will contribute to our knowledge of the underlying 
causes of ASD-related behaviour. Previous research high- 
lights the importance of investigating cognitive impair- 
ments in ASD that may arise from brain abnormalities 
and drive behavioural symptoms [54]. Further, exploring 
WM correlates will supplement behavioural phenotypes of 
ASD. This work will help identify the nature of atypical 
brain development, with future expectations of establish- 
ing age-appropriate interventions that can effectively tar- 
get WM function and, in turn, other symptoms of ASD. 

Methods 

Participants 

Seventy-three participants were recruited for this study: 
42 children aged 7-13 years with high-functioning ASD, 
and 31 7-13-year-old TD control children. However, 9 
TD children and 23 children with ASD were excluded 
from the analyses for excessive movement and inadequate 
task performance or protocol completion. After sex-, 
IQ- and age-matching, the study sample consisted of 19 
children with ASD (3 girls and 16 boys) and 17 controls 
(4 girls and 13 boys). The groups were matched for age 
(ASD M= 11.05, SD = 1.43; TD M = 11.12, SD = 2.00; 
£(34) = 0.11, ns.), sex (xfi) = 0.34, ns.) and full-scale IQ 
as determined by the Wechsler Abbreviated Scale of 
Intelligence-II [55] (ASD M = 109.42, SD = 15.72; TD 
M = 115.35, SD = 9.27; t m = 1.36, ns.). We substituted the 
group average IQ for one TD child whose data were 
missing. 

Exclusion criteria for all participants were the presence 
of any current significant axis I psychiatric comorbidities 
[2], neurological disorders, medical illnesses, prematurity, 
uncorrected vision, colour blindness, IQ < 70 as well as 
standard MRI contraindicators (e.g. ferromagnetic im- 
plants). A history of developmental delay, learning disabil- 
ity and attention deficit hyperactivity disorder (ADHD) 
was used to exclude TD children only; however, these 
factors were also not current primary diagnoses in the 
ASD group. Six children with ASD were each on one 
psychotropic medication (Strattera, Biphentin, Fluoxetine, 
Concerta, Ability and Atomoxetine). Their fMRI data were 
examined in comparison to children with ASD who were 
not taking medication, and the data did not differ between 
these subgroups (see Additional file 1). Children were 
recruited through community support centres, parent 
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support groups, email lists, hospital ads and private 
schools. Informed consent, clinical and cognitive testing 
and MRI scanning were performed at the Hospital for 
Sick Children in Toronto. Experimental procedures were 
approved by the Research Ethics Board at the Hospital for 
Sick Children. All children gave informed assent, and the 
parents provided informed written consent. 

Clinical diagnoses of ASD were confirmed in all cases 
with a combination of expert clinical judgement, clinical 
records and the Autism Diagnostic Observation Schedule 
(ADOS) [56], which was administered by a trained indi- 
vidual who maintains inter-rater research reliability. All 
children completed the backwards digit recall, listening 
recall, digit recall, mazes memory and block recall sub- 
scales of the Working Memory Test Battery for Children 
(WMTB-C) [57] to supplement behavioural data collected 
during fMRI tasks. See Table 1 for demographic, neuro- 
psychological and clinical test characteristics. 

The colour matching task 

Children were required to attend to coloured figures of a 
clown presented one at a time in sequence. Children 
were taught to ignore the clown's face and irrelevant 
colours (blue and green) and focus on other, relevant 
colours (yellow, purple, pink, red, orange, brown and 
grey). Tasks that contain misleading or irrelevant factors 
that evoke interference have been shown to be a more 
suitable measure of WM capacity, yielding reliable esti- 
mates of developmental growth in adults [48] and children 
(Powell, Arsalidou, Vogan, and Taylor: Controlled infer- 
ence and assessments of developmental working memory 
capacity: Evidence from letter and colour matching tasks, 
submitted). The number of 'n relevant colours (capacity) 
in the figures was increased by one for each increase in 
difficulty level. CMT has two integral functions that re- 
quire mental attention, in which participants must first 



actively extract the relevant colours embedded in the 
clown figure and second check for a possible match 
with colours of the criterion set. As such, items with n 
(e.g. 1) relevant colours will have difficulty level of n + 2 
(e.g. 3; Figure 1A). Children indicated after each item 
whether the relevant colours of the current figure 
matched those from the immediately preceding figure 
(i.e. 1-back), disregarding colour location and repetition. 
Using a keypad with the right hand, children responded 
in the scanner by pushing a button for 'same' when the 
clown figure was wearing the same colours as the previ- 
ous clown and 'different' when the clown figure was 
wearing different colours. All children successfully com- 
pleted practice trials on a computer outside the scanner 
with accuracy of 80% or greater. 

A total of 24 task blocks (168 task trials) and 24 base- 
line blocks were presented across four runs. Each run 
consisted of six 32-s blocks, one for each difficulty level 
(six levels in total); within each block, there was a constant 
difficulty level, and all difficulty levels were randomized 
within each run (Figure IB). The same four runs were pre- 
sented to all children in the same order. Task blocks con- 
sisted of eight stimuli of the same difficulty and alternated 
with 20-s baseline blocks (Figure 1C), where clowns were 
coloured only in blue and green, and children were 
instructed to look at the figures but not respond. Partici- 
pants had 3 s to view a stimulus and respond, followed by 
a 1-s inter-stimulus interval where a fixation cross was 
presented (Figure ID). The fMRI task took approximately 
22 min of scan time. 

Performance data were recorded for both accuracy 
(proportion correct) and reaction time; items were correct 
if responded correctly within 3 s of stimulus presentation. 
To ensure adequate task completion, children were ex- 
cluded if they did not achieve at least 60% accuracy (aver- 
aged across four runs) on the easiest two difficulty levels, 



Table 1 Demographic and neuropsychological test characteristics of study sample 

Variables ASD(A/=19) TD(N=17) 



% 



Mean (SD) 



% 



Mean (SD) 



Significance test 



Demographic data 
Sex (male) 
Age 

IQ 

ADOS Total 
Neuropsychological test data (WMTB-C) 
Digit recall 
Block recall 
Mazes memory 
Listening recall 
Backward digit recal 



84.21 



11.05 (1.43) 
109.42 (15.72) 
9.68 (2.3) a 

107.59 (19.53) 
95.58 (18.03) 
95.42 (18.58) 
102.95 (15.77) 
94.05 (18.79) 



76.47 



11.12 (2.00) 
115.35 (9.27) 
N/A 

1 1 6.00 (1 5.64) 
105.35 (20.72) 
96.35 (2349) 
116.76 (17.72) 
109.18 (21.09) 



f m = 0.34, p = 0.56 
f (34) = 0.1 1,p = 0.91 
r (34) = 1.36, p = 0.18 



f <34) = 1 .51, p = 0.14 
t (34) = 1.51, p = 0.14 
f (34) = 0.13, p = 0.90 
r (34) = 2.48, p = 0.02* 
f (34 , = 2.28, p = 0.03* 



*p < 0.05. a ADOS scores can range from 3-20, with higher scores reflecting greater symptom severity. 
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Figure 1 (See legend on next page.) 
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(See figure on previous page.) 

Figure 1 Protocol description of the colour matching task (CMT). (A) There were six levels of difficulty where the number of relevant colours 
(yellow, purple, pink, red, orange, brown and grey) increased by one to increase the difficulty level. Difficulty = (# of colours) + 2. Children were 
taught to ignore the clown's face, colour location, colour repetition and irrelevant colours (blue and green). (B) It was a block design task, where 
each run consisted of six 32-s task blocks (for each difficulty) followed by 20-s baseline blocks where clowns are presented in only blue and green 
(ignore). Task blocks were presented pseudo-randomly within each run. (C) Example of part of a sequence in a baseline block. Stimuli were 
presented for 3 s followed by 1-s inter-stimulus cross. Children were instructed to not respond. (D) Example of part of a sequence in a task block; 
children indicated if the current clown was wearing the same or different colours as the previous clown. Stimuli were presented for 3 s followed 
by a 1-s inter-stimulus fixation cross. 



and also excluded if they did not have at least two out of 
four runs where 50% or more of the blocks were accept- 
able in terms of performance (60% accuracy) and motion. 
Motion was deemed acceptable if children moved less 
than 1.5 mm from their median head position in at least 
60% of the volumes within a task block. See the fMRI pre- 
processing section below for a description of displacement 
calculations. A 60% accuracy criterion was chosen because 
while we could be sure that children were performing bet- 
ter than chance (50%), it was not too stringent. Motion pa- 
rameters were also entered into the fMRI preprocessing 
pipeline. 

Image acquisition 

All imaging data were acquired using a 3 T Siemens 
Trio MRI scanner with a 12-channel head coil. Head 
stabilization and motion restriction were achieved with 
foam padding. The structural scan was a high-resolution 
Tl-weighted 3D MP-RAGE image (Sagittal; FOV = 192 x 
240 x 256 mm; 1 mm isometric voxels; TR/TE/TI/FA = 
2,300/2.96/900/9), which was used as an individual ana- 
tomical reference for the functional images. During 
structural image acquisition, children watched a movie 
of their choice using MR-compatible goggles and ear- 
phones. Functional images were acquired with single-shot 
echo planar imaging sequence (Axial; FOV = 192 x 192; 
Res = 64 x 64; 30 slices 5 mm thick; 3x3x5 mm voxels; 
TR/TE/FA = 2,000/30/70). Visual stimuli for the functional 
task (CMT) were displayed on MR-compatible goggles. 
Children responded to trials using a dual button MR- 
compatible keypad. Stimuli were displayed and per- 
formance was recorded using the software Presentation 
(Neurobehavioral Systems Inc., Berkeley, CA, USA). 

Behavioural data analyses 
CMT 

Accuracy on difficulty levels 7 and 8 (D7 and D8) was 
poor for both TD children (D7, M = 0.57, SD = 0.03; D8, 
M = 0.53, SD = 0.02) and children with ASD (D7, M = 
0.53, SD = 0.02; D8, M = 0.54, SD = 0.03); thus analyses 
of only the first four difficulty levels (D3 to D6) were 
completed. Accuracy and response times were calculated 
for each difficulty level by averaging across runs for each 
group. Data were analyzed using repeated measures 



factorial ANOVAs, with group (ASD and TD) as a be- 
tween subject factor and difficulty level (D3, D4, D5 and 
D6) as a within subject factor. 

WMTB-C 

Standardized scores on the subscales were compared 
across group using t tests to determine if there were dif- 
ferences between ASD and TD children on these neuro- 
psychological measures of WM. 

fMRI data analyses 

Image preprocessing of functional data was performed 
using a combination of standard AFNI [58] and FMRIB's 
Software Library (FSL) [59] tools. The first three volumes 
of each run were discarded for scanner stabilization. After 
slice timing and motion correction, data were smoothed 
using a 6-mm FWHM Gaussian kernel, temporally filtered 
(lower and upper cutoff frequencies of 0.01 and 0.2 Hz, re- 
spectively) and converted to percent signal change from 
the baseline volumes. Before group-level analyses, images 
were registered to the Montreal Neurological Institute 
(MNI) 152 template. The maximum Euclidean displace- 
ment (MD) travelled by any voxel within the brain was 
calculated from the six rigid body transformation parame- 
ters for each volume. This MD metric was used to flag vol- 
umes with unacceptable motion, as described above. The 
average MD for each subject was used to explore group dif- 
ferences in head motion. Although more motion was found 
in children with ASD (M = 0.58 mm, SD = 0.54 mm) than 
TD children (M = 0.29 mm, SD = 0.22 mm), t {34) = 2.07, 
p = 0.05, both groups had minimal average motion of 
under 0.60 mm. To control for motion, MD was also in- 
cluded as a covariate of no interest in the GLM. 

Data were analyzed with the FSL fMRI Expert Analysis 
Tool (FEAT) [60]. Data were fit first to a block-design 
general linear model convolved with a gamma function 
to model haemodynamic response, using the task parame- 
ters (D3 to D6). To examine areas that linearly modulated 
as a function of difficulty, linear trend analyses were con- 
ducted from D3 to D6 using fixed-effects higher level 
modelling. Individual results were then averaged across 
runs for each subject in a second-level analysis. Between- 
group comparisons were carried out using FMRIB's Local 
Analysis of Mixed Effects- 1 (FLAME 1) [59] to obtain an 
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accurate between-subject variance estimation, which in- 
creased our ability to detect real activation [60] . Significant 
activations were reported using cluster-based thresholding 
determined by Z > |2.3| and a corrected cluster signifi- 
cance threshold of p < 0.05. Regions of interest (ROIs) 
were selected from the local maxima of areas showing sig- 
nificant group differences between TD and ASD groups in 
the linear trend analyses. Average percent signal change 
and standard error scores were extracted from spherical 
ROIs (6 mm radius) centred about the local maxima of 
group difference maps, and the mean peak cluster signal 
change for each group was plotted as a function of diffi- 
culty to further examine activation patterns across cogni- 
tive load. 

Results 

Behavioural data 

There was a significant effect of group on accuracy 
(Figure 2A), F (1 , 34) = 5.15, p = 0.03, n 2 = 0.13, which was 
driven by TD children (M = 0.80, SD = 0.12) performing 
more accurately than ASD children (M = .71, SD = 0.12) 
only at D5 (t(34) = 2.30, p = 0.03). As such, comparisons of 
brain activity between control and ASD groups were made 
under comparable accuracy scores across most levels. 
There was a significant main effect of difficulty on accur- 
acy, F(3, 34) = 53.87, p < 0.001, tf = 0.61, with performance 
accuracy decreasing as a function of difficulty level in both 
groups. Post hoc pairwise comparisons, adjusted for mul- 
tiple comparisons using Bonferroni, revealed that accuracy 
on the most difficult level was significandy different from 
each other at p < 0.05, the exceptions being between D5 
and D6 in the ASD group (Table 2), and between D3 and 
D4, and between D4 and D5 in TD children (Table 2). We 
also ran a supplementary analysis, inputting performance 
as a covariate. There were no areas of activation that cor- 
related with performance in either group, suggesting that 



findings were not confounded by performance per se. 
Also, there were no significant differences in performance 
across runs in the TD (i r ( 3j 48) =1.21, p = 0.32) or ASD 
(F (3, 45) = 0.50, p = 0.69) group, and therefore patterns 
of performance did not reflect fatigue across time. 

Overall, response times increased with increasing diffi- 
culty in TD children but only increased up until D5 in 
children with ASD (Figure 2B). There was no main effect 
of group on response times (F (lj 34) = 2.33, p = 0.14), but 
there was a main effect of difficulty level, F( 3i 34) = 54.60, 
p < 0.001, rj 1 = 0.62. Post hoc comparisons in the ASD 
group showed that response times differed between dif- 
ficulty levels, except between D4 and D6, and between 
D5 and D6 (Table 3). In TD children, response times 
significantly differed between difficulty levels, except 
between D4 and D6 and D5 and D6 (Table 3). 

For the WMTB-C, TD children had significantly higher 
scores on the listening recall (f( 34 ) = 2.48, p < 0.05) and 
backward digit recall (i( 34 ) = 2.28, p < 0.05) subtests than 
children with ASD. The groups did not differ on any other 
subtests (see Table 1 for scores on the WMTB-C). 

fMRI Data 

Task-related activation within groups 

The primary objectives of the fMRI analyses were to inves- 
tigate the pattern of brain activity exhibited as a function 
of cognitive load (i.e. difficulty level) and determine if this 
pattern differed in children with and without ASD. Linear 
trend analyses (D3 to D6) showed that while some brain 
areas increased in activity as a function of difficulty level, 
others decreased. 'Increasing activation' refers to an in- 
crease in BOLD signal with increasing load (i.e. positive 
linear relations between cortical activity and task difficulty 
level) and 'decreasing activation' refers to a decrease in 
BOLD signal with increasing load (i.e. negative linear rela- 
tions between cortical activity and task difficulty level). As 
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Figure 2 CMT behavioural performance. (A) Mean proportion correct for D3 to D6 and standard error bars. TD children were significantly 
more accurate than children with ASD at D5 only. (B) Mean response times for D3 to D6 and standard error bars. There were no significant 
differences between TD and ASD groups in response times across difficulty levels. 
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Table 2 Differences in CMT accuracy (proportion correct) 
across difficulty 



Differences in CMT accuracy 



D3 




D4 




D5 




MD 


SE 


MD 


SE 


MD 


SE 


Children with ASD 












D4 0.06** 


0.02 










D5 0.17*** 


0.03 


0.11** 


0.03 






D6 0.23*** 


0.03 


0.17*** 


0.03 


0.07 


0.03 


TD Children 












D4 0.05 


0.02 










D5 0.11* 


0.03 


0.07 


0.02 






06 0.22*** 


0.03 


0.17*** 


0.03 


0.11** 


0.03 



*p < 0.05, **p < 0.01, ***p < 0.001, corrected for multiple comparisons using 
Bonferroni. Posf hoc ANOVA tests using accuracy scores for (A) children with 
ASD and (B) TD children. D3 to D6 difficulty levels 3 to 6, MD mean difference, 
SE standard error. 



shown in Figure 3, magnitude of the signal change in- 
creased with difficulty in TD children, particularly be- 
tween D5 and D6, in bilateral fusiform (BA37), precuneus 
(BA7), inferior frontal gyri (BA45), right DLPFC (BA9) 
and bilateral dorsal cingulate/dorsal medial prefrontal cor- 
tex (B32/8) extending to the anterior cingulate (B33/24). 
In the bilateral middle occipital gyrus (BA19), a positive 
linear change in activation with increasing task difficulty 
was seen up to D6 (Table 4). Activation found in the bilat- 
eral posterior cingulate (BA 23/31) and anterior medial 
prefrontal gyrus (BA10) decreased (i.e. showed a negative 
linear relation) as a function of difficulty level (Table 4). 

Children with ASD did not show the same trends in 
activation, particularly in the frontal brain regions. A 
positive linear relation between cortical activity and task 
difficulty was seen only in the posterior brain regions, 

Table 3 Differences in CMT response times (seconds) 
across difficulty levels 



Differences in CMT response times 





D3 




D4 




D5 




MD 




SE 


MD 


SE 


MD 


SE 


Children with ASD 












D4 -0.25* 




0.05 










D5 -0.38* 




0.05 


-0.13** 


0.03 






D6 -0.33* 




0.07 


-0.08 


0.07 


0.05 


0.05 


TD Children 














D4 -0.22* 




0.03 










D5 -0.41* 




0.04 


-0.19*** 


0.03 






D6 -0.47* 




0.07 


-0.25** 


0.06 


-0.06 


0.05 



*p < 0.05, ** p < 0.01, ***p < 0.001, corrected for multiple comparisons using 
Bonferonni. Post hoc ANOVA tests using response times for (A) children with 
ASD and (B) TD children. D3 to D6 difficulty levels 3 to 6, MD mean difference, 
SE, standard error. 



including the bilateral middle occipital (BA19) and fusi- 
form gyri (BA37, Table 5). Similar to TD children, acti- 
vation in the posterior cingulate (BA23/31) and medial 
prefrontal gyrus (BA 10/32) decreased (i.e., showed a 
negative linear relation) as a function of difficulty level 
in children with ASD (Table 5). 

Between-group differences in task-related activation 

Between-group analyses of each difficulty level separately 
showed no significant differences (see Additional file 2). 
However, there were significant group effects in linear 
activation trends as a function of difficulty level. As 
shown in Figure 4, three regions demonstrated group 
differences in activation. These areas included the bilat- 
eral precuneus (BA7), right DLPFC (BA9) and left dorsal 
medial premotor cortex (BA8) (Table 6). In these regions, 
TD children showed significant positive linear relations 
between cortical activity and difficulty level, particularly 
between D5 and D6, whereas children with ASD did not 
show these trends between D5 and D6. The magnitude of 
the signal change in frontal and parietal areas increased 
largely with the most difficult level when TD children 
tried to meet WM demand, but children with ASD failed 
to show this increase. See Figure 5 for graphs of percent 
signal change with standard error for ROIs of brain re- 
gions that showed significantly different linear patterns 
between TD and ASD groups. 

Discussion 

This is the first study to examine the neural correlates of 
visuo-spatial WM in pre-adolescent children with ASD 
relative to TD children. Whereas children with and with- 
out ASD did not differ in brain activation during WM 
function irrespective of load, between-group differences 
were observed in the linear activation trends across diffi- 
culty level. In other words, children with and without 
ASD differ in how they modulate WM processes during 
tasks that increase in difficulty level. Using a task that iso- 
lated cognitive load on WM, we observed positive linear 
relations between cortical activity and task difficulty level 
in prefrontal and parietal regions, particularly between the 
highest difficulty levels, in TD children that were signifi- 
cantly different from activation trends seen in children 
with ASD. These areas included bilateral precuneus, right 
DLPFC and left medial premotor cortices. Other areas 
known to underlie WM function, including the anterior 
cingulate and inferior frontal gyri, followed this same pat- 
tern but the linear activation trends did not differ signifi- 
cantly between TD and ASD subjects, likely due to subject 
variability. Overall, TD children demonstrated an oppos- 
ing system of cognitive processes where areas related to 
task difficulty (frontal regions) increased in activity and 
areas associated with the brain's DMN (posterior cingulate 
and anterior medial frontal gyrus) decreased in activity 
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Figure 3 Group activation maps for the linear trend analyses in ASD and TD groups during CMT. Significant activations using cluster-based 
thresholding determined by Z > |2.3| and a corrected cluster significance threshold of p = 0.05. Areas in red depict regions of increasing activation as a 
function of difficulty (i.e. positive linear relations between cortical activity and task difficulty level), and areas in blue depict regions of decreasing 
activation (i.e. negative linear relations between cortical activity and task difficulty level). dmPFC dorsal medial prefrontal cortex, Post Cing 
posterior cingulate cortex, Prec precuneus, DLPFC dorsolateral prefrontal cortex, ACC anterior cingulate cortex, IFG inferior frontal gyrus, MidOcG 
middle occipital gyrus, FusG fusiform gyrus, amPFG anterior medial prefrontal gyrus. 
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Table 4 Linear trend analyses across difficulty levels: TD group 



Linear trend analyses across difficulty levels for TD children 





Voxels 


MNI Coordinates 


Z value 


P value 


Hem. 


Region 






X 


y 


z 










Regions where activation increases with 


8885 


-26 


-90 


16 


5.04 


1.88 X 10~ 21 


L 


Middle occipital gyrus 


difficulty (increasing BOLD signal) 


X 


24 


-82 


-10 


4.93 




R 


Fusiform/lingual gyrus 




X 


24 


-62 


50 


4.86 




R 


Precuneus 




V 
A 


38 


-84 


10 


A 19. 

H. 1 O 




D 

r\ 


• • 

Middle occipita gyrus 




X 


-16 


-64 


56 


4.25 




L 


Precuneus 




X 


-20 


-82 


-14 


4.22 




L 


Fusiform/lingual gyrus 




6845 


32 


26 


0 


4.97 


7.26 X 10~ 18 


R 


Inferior frontal gyrus 




X 


-32 


24 


-4 


4.95 




L 


Inferior frontal gyrus 




X 


8 


28 


32 


4.94 




R 


Anterior cingulate cortex 




X 


32 


42 


22 


4.04 




R 


Dorsolateral prefrontal cortex 




X 


-12 


30 


24 


3.90 




L 


Anterior cingulate cortex 


Regions where activation decreases with 


1808 


-6 


-56 


30 


-3.68 


9.54 X 10"' 


L 


Posterior cingulate gyrus 


difficulty (decreasing BOLD signal) 


X 


4 


-40 


34 


-3.47 




R 


Posterior cingulate cortex 




725 


-42 


-72 


36 


-3.69 


3.29 X 1 0" 3 


L 


Middle temporal gyrus 




629 


8 


'18 


42 


-3.76 


7.91 X 10 3 


R 


Medial frontal gyrus 




X 


-6 


62 


20 


-3.61 




L 


Medial frontal gyrus 



Results from linear trend analyses from D3 to D6 for TD children. Areas that increased as a function of difficulty level (A) are associated with WM and visuo-spatial 
processing, whereas areas that decreased as a function of difficulty level (B) are associated with the default mode network. MNI coordinates represent the peak Z 
value of the cluster; X peak local maximas within cluster. 



with increasing cognitive demand; this pattern of activa- 
tion was absent in children with ASD. There were signifi- 
cant differences in performance on only two verbal WM 
subtests of the WMTB-C. Further, TD and ASD groups 
performed similarly on the CMT (except for a small differ- 
ence at D5), suggesting that the findings are not con- 
founded by behavioural differences. 

TD children showed increased recruitment of both 
frontal and posterior parietal regions as a function of 
cognitive load, whereas the ASD group showed in- 
creased recruitment of only the posterior visual regions, 
including the left and right fusiform and middle occipi- 
tal gyri. Our observed pattern of findings in children 
with ASD are consistent with findings in adults with 
ASD who also relied on posterior regions related to 
low-level cognition during WM tasks, rather than areas 
of high level-cognition, such as in the prefrontal corti- 
ces [10]. More activation in occipital-temporal areas in 
ASD has been proposed to reflect visually based pro- 
cessing styles and a tendency to rely primarily on visual 
features and details of objects, rather than on WM [61]. 
These findings are in line with the enhanced perceptual 
functioning model [62] that suggests individuals with 
ASD display superior activation of visuo-perceptual re- 
gions in association with a reduced activation in higher 
order frontal areas. The present study extends this lit- 
erature in adults by highlighting that this processing 
style appears in young children with ASD as well. The 



same pattern of activation has also been found in so- 
cially relevant tasks [63,64]. Although greater posterior 
activation may be adequate for visuo-spatial processing, 
it is likely inefficient for more complex cognitive opera- 
tions, such as language comprehension [19] and WM 
[10], which would have implications for the social and 
executive deficits typically observed in ASD. Less exten- 
sive use of prefrontal areas may be the result of early ab- 
normal prefrontal development [15,16] and connectivity 
[10,19] reported in individuals with ASD. 

The finding of stable activation with increasing cognitive 
load within the prefrontal cortex may be most significant 
for ASD due to the crucial role of frontal regions in WM 
and executive function [10,43,44,65-67]. Similar to adults 
[10,12] and adolescents [14] with ASD, we found that 
pre-adolescent children with ASD demonstrated differ- 
ent activation patterns in the DLPFC (BA 9) during 
WM function. Specifically, in the DLPFC, TD children 
showed a graded increase in activity from D3 to D6, 
whereas activity within this region did not increase at 
the same magnitude in children with ASD. The DLPFC 
is believed to play a critical role in holding information 
'online' [26] and mediating strategic organization and 
data compression processes [68,69], hence its sensitivity 
to increasing cognitive demand in ours and other stud- 
ies [68-71]. Stable activation across difficulty level in 
this area suggests that young children with ASD may 
fail to use appropriate organizational strategies, such as 
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Table 5 Linear trend analyses across difficulty levels: ASD group 



Linear trend analyses across difficulty levels for children with ASD 





Voxels 


MNI Coordinates 


Z value 


P value 


Hem. 


Region 






X 


y 


z 










Regions where activation increases with 
difficulty (increasing BOLD signal) 


1816 


28 


-52 


-14 


3.61 


8.94 X 1 0~ 7 


R 


Fusiform gyrus 




v 

A 


26 


-84 


12 


3.42 




n 

n 


Midd e occipital gyrus 




1449 


-26 


-88 


8 


3 52 


1 09 X 1 0 -5 


|_ 


(\/liHH p r™ii~i~i nita n\/n 
IVIIUUIt: Ul_l_IUILGI yyiui 




X 


-20 


-68 


-14 


3.17 




L 


Fusiform gyrus 


Regions where activation decreases with 
difficulty (decreasing BOLD signal) 


2210 


10 


-50 


36 


-3.80 


5.96 x 10~ 8 


R 


Posterior cingulate cortex 




X 


-6 


-68 


24 


-2.91 




L 


Posterior cingulate cortex 




1308 


12 


54 


0 


-4.06 


3.02 X 10~ 5 


R 


Medial frontal gyrus 




X 


-10 


40 


8 


-2.98 




R 


Medial frontal gyrus/Anterior cingulate cortex 




642 


18 


36 


44 


-3.64 


7.00 X 10" 3 


R 


Superior frontal gyrus 




509 


60 


-56 


-2 


-3.63 


2.89 X 10~ 3 


R 


Middle temporal gyrus 



Results from linear trend analyses from D3 to D6 for children with ASD. Areas that increased as a function of difficulty level (A) are associated with WM and 
visuo-spatial processing, whereas areas that decreased as a function of difficulty level (B) are associated with the default mode network. MNI coordinates represent 
the peak Z value of the cluster; X peak local maximas within cluster. 



'chunking' methods, which facilitate WM by simplifying 
cognitive load/demand. Furthermore, normative studies 
have demonstrated that with increasing age, individuals 
rely more on DLPFC in WM processes [43,66,72], sug- 
gesting that specialization of this region for WM coin- 
cides with structural maturation across development. It 
is possible that atypical abnormal growth patterns in 
ASD in the frontal cortex, the DLPFC in particular [15], 
adversely impacts its functional integrity; these specula- 
tions require further research. 



Failure to increase recruitment of the precuneus across 
difficulty level in ASD subjects is also of interest, given 
evidence supporting parietal involvement in addition 
to the prefrontal cortex during spatial WM processing 
[44,66,70,73]. The precuneus forms part of the proposed 
occipito-parietal network, or alternatively the visual 'dorsal 
pathway', responsible for spatial visual processing (e.g. 
object location) [74]. Conversely, the occipitotemporal, 
or visual 'ventral pathway! includes the fusiform gyri and 
is critical to object identification (e.g. colour and shape) 



Prec (BA 7) 



DLPFC 
(BA 9) 



dmPMC 
(BA 8) 



Prec 
(BA 7) 



Figure 4 Results from between-group comparisons. Significant activations using cluster-based thresholding determined by Z> 2.3 and a 

corrected cluster significance threshold of p = 0.05. Areas in red depict regions where the control children showed greater positive linear relations 
between cortical activity and task difficulty level than the ASD group. There were no areas where the ASD group showed greater linear activation 
trends across difficulty level in the negative or positive direction than controls. Prec precuneus, DtPfC dorsolateral prefrontal cortex, dmPMC dorsal 
medial premotor cortex. 
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Table 6 Regions of significant differences between TD and ASD groups 


Voxels 




MNI Coordinates 




Z value 


P value 


Hem. 


Region 




X 


y 


z 










810 


0 


-56 


46 


3.93 


1.56X 10~ 3 


L 


Precuneus 


X 


10 


-66 


40 


3.81 




R 


Precuneus 


796 


34 


34 


40 


3.65 


1.76X 10~ 3 


R 


Dorsolateral prefrontal cortex 


X 


-6 


26 


46 


3.13 




L 


Dorsal medial premotor cortex 



Results from between group comparisons of the linear trend analyses from D3 to D6. All regions reported are areas in which TD children showed greater positive 
linear relations between cortical activity and difficulty level (i.e. increasing BOLD signal with increasing difficulty) than children with ASD. There were no areas 
where children with ASD showed greater linear relations between cortical activity and difficulty level than TD children. MNI coordinates represent the peak Z 
value of the cluster; X peak local maximas within cluster. 



[74]. Given that CMT was a visuo-spatial task with 
gradual increases in WM that needed to be processed, 
areas responsive to both spatial search (occipito-parietal 
regions) and categorization (occipito-temporal regions) 
should show increased recruitment across cognitive load. 
Although group differences in fusiform activity were 
absent, we observed a significant difference in the pre- 
cuneus; TD children showed positive linear relations 
between cortical activity and difficulty level, particu- 
larly between D5 and D6, whereas the ASD group did 
not show this linear change. Impaired parietal activity 
in pre-adolescent children with ASD lends support to 
the proposed dorsal stream deficits within this population 
as well as in other developmental disorders [75-78], 



while ventral processing is relatively intact. Further- 
more, normative developmental studies suggest that the 
dorsal pathway has a more protracted maturational tra- 
jectory than the ventral stream [66], and WM fMRI 
studies demonstrate greater parietal [4,73] but less fusi- 
form recruitment [66] with age. Thus, spared fusiform/ 
occipital activity and abnormal parietal functioning in 
our ASD group may reflect immature WM processing 
typically seen in very young TD children. However, it is 
not clear whether the recruitment of mature neural sub- 
strates underlying WM processes is delayed, persistently 
weak or arrested. Future longitudinal work will help 
clarify the developmental path of WM neural circuitry 
in ASD. 
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Figure 5 Mean peak cluster percent signal changes and standard error. As a function of difficulty between task difficulty and baseline conditions 
in areas where children with ASD significantly differed from TD children in the linear trend analyses. DLPFC dorsolateral prefrontal cortex. 
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Conclusions 

Overall, the current study fills a significant gap in our 
knowledge of neural substrates implicated in visuo-spatial 
WM functioning among pre-adolescent children and how 
they differ in children with ASD. Similar to adults, we 
found that children with ASD show prefrontal and parietal 
system abnormalities and tend to rely on posterior brain 
regions associated with lower level cognitive processing. 
The present study also converges on the growing body of 
literature proposing dorsal stream deficits in visuo-spatial 
processing in this population. This study contributes new 
information on WM differences in school-aged children 
with and without ASD, demonstrating that even in pre- 
adolescents, there are significant differences in brain acti- 
vation patterns during WM processing with increasing 
cognitive load that differentiate the groups. In light of 
similar task performance between groups, it is also im- 
portant to entertain the idea that higher order processing 
may not be mandatory in ASD with tasks that can be 
processed efficiently using a perceptual processing ap- 
proach because of their proposed enhanced perceptual 
functioning [62] . Conversely, higher order control may be 
necessary in TD children who do not have such perceptual 
expertise. Although data may support this argument, ASD 
is generally considered a network/brain disorder, and what 
is interesting is how networks in children with ASD differ 
from TD children with complex tasks. Group effects were 
driven by the most difficult level analyzed, which may be 
due to strategies used by TD children necessary to encode 
the increased number of colours in the highest task load 
(i.e. meet high WM demand) that children with ASD are 
unable to employ. Similarly, previous research reports be- 
havioural differences with tasks that have more complex 
WM demands e.g. [9,34,39]. 

It is important to consider the limitations of the current 
study when interpreting results. With comparable CMT 
behavioural performance between ASD and TD groups, 
we eliminated performance as a confounding factor. Con- 
sequently, our sample was less representative of low- 
functioning individuals with ASD, and thus, results are 
generalizable to higher functioning individuals only. Fu- 
ture fMRI research is required to understand WM func- 
tion across various levels of functioning and a range of 
symptoms. In addition, given our choice of control sub- 
jects (TD children), findings can only provide informa- 
tion about differences from the norm. Comparisons to 
other atypical populations who share similar cognitive 
but different clinical profiles (e.g. ADHD) [5] will further 
our understanding about the neural patterns that are 
unique to ASD, potentially explaining characteristic behav- 
iour in this complex group. Lastly, results are also limited 
by the relatively small sample size, due to assessing a com- 
plex cognitive ability while scanning, which greatly in- 
creased the amount of children's data that was not usable. 



The findings from the current study have a number of 
significant implications despite these limitations. Several 
researchers have stressed the importance of executive 
cognitive skills for social function [5,79], a core deficit of 
ASD. However, there remains a gap in knowledge regard- 
ing the link between neuropsychopathology and clinical 
symptoms of autism. With respect to our findings of atyp- 
ical neural activity underlying WM processes in ASD, this 
may impair the ability to hold information 'online' that 
may affect one's ability to evaluate and select appropriate 
responses during peer interactions, translating into socially 
inappropriate behaviour. Social behaviour requires com- 
plex cognitive processing, further highlighting the import- 
ance of our results of increasing group differences in brain 
activity with increases in complexity. Future neuroimaging 
studies could work towards understanding the relation be- 
tween ASD symptomology and neural activation patterns 
associated with WM processing. Overall, our findings will 
help guide future longitudinal work by localizing areas of 
vulnerability to developmental disturbances and allow 
health care providers to carefully monitor their develop- 
ment. Developmental information will allow us to iden- 
tify the nature and timing of atypical development, 
which is critical in establishing age-appropriate cogni- 
tive or pharmacological remediation for WM function 
and behavioural deficits in ASD. 

Additional files 



Additional file 1: fMRI data for medicated and non-medicated 
children with ASD. Percent signal change as a function of difficulty 
between task difficulty and baseline conditions in children with ASD who 
were on medication versus those who were not. Areas of the brain 
shown are from regions where children with ASD significantly differed 
from TD children in the linear trend analyses. Potential differences 
between medicated and non-medicated children with ASD were also 
examined statistically using the FSL FEAT, and no significant differences 
were found between children with ASD who were and were not on 
medication. However, due to low N (only six subjects on medication), this 
statistical test may not be reliable. Therefore scatter plots were created to 
visually examine the data for significant group differences, and this 
reaffirmed that medication does not appear to affect the findings. 

Additional file 2: Individual group activation maps for children with 
ASD and typically developing children at all levels (D3 to D6) of 
CMT. Significant activations using cluster-based thresholding determined 
by Z> 1 2.3 1 and a corrected cluster significance threshold of p = 0.05. 
Areas in red and blue depict regions with significantly higher and lower 
BOLD signal than baseline, respectively. Between-group comparisons 
showed no areas of significant difference between children with and 
without ASD at any single difficulty level. 
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